Projected Constraints on Modified Gravity Cosmologies 
from 21 cm Intensity Mapping 
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We present projected constraints on modified gravity models from the observational technique 
known as 21 cm intensity mapping, where cosmic structure is detected without resolving individual 
galaxies. The resulting map is sensitive to both BAO and weak lensing, two of the most powerful 
cosmological probes. It is found that a 200 m x 200 m cylindrical telescope, sensitive out to z = 2.5, 
would be able to distinguish DGP from most dark energy models, and constrain the Hu & Sawicki 
f{R) model to |/ho| < 9 x 10 -6 at 95% confidence. The latter constraint makes extensive use of 
the lensing spectrum in the nonlinear regime. These results show that 21 cm intensity mapping is 
not only sensitive to modifications of the standard model's expansion history, but also to structure 
growth. This makes intensity mapping a powerful and economical technique, achievable on much 
shorter time scales than optical experiments that would probe the same era. 
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I. INTRODUCTION 

One of the greatest open questions in cosmology is the 
cause of the observed late time acceleration of the uni- 
verse. Within the context of normal gravity described by 
Einstein's General Relativity, this phenomena can only 
be explained by an exotic form of matter with negative 
pressure. Another possible explanation is that on cos- 
mological scales, General Relativity fails and must be 
replaced by some theory of modified gravity. 

Several approaches have been proposed to modify grav- 
ity at late times to explain the apparent acceleration of 
the universe. The challenge in these modifications is to 
preserve successful predictions of the CMB at z 1000, 
and also the precision tests at the present epoch in the 
solar system. 

A generic class of theories operates with the 
Chameleon effect, where at sufficiently high densities 
General Relativity (GR) is restored, thus applying both 
in the solar system and the early universe. To further 
understand the nature of gravity would require prob- 
ing gravity on cosmological scales. Large scales means 
large volume, requiring large fractions of the sky. Gravity 
can be probed by gravitational lensing, which measures 
geodesies and thus the gravitational curvature of space, 
and is a sensitive probe of the growth of structure in the 
Universe 

In working out predictions for cosmology, the theo- 
retical challenge posed by these theories are the non- 
linear mechanisms in each model, necessary in order to 
restore Einstein Gravity locally to satisfy Solar System 



constraints. We present quantitative results from nonlin- 
ear calculations for a specific f(R) model, and forecasted 
constraints for future 21cm experiments. 

An upcoming class of experiments propose the obser- 
vation of the 21 cm spectral line at low resolution over 
a large fraction of the sky and large range of redshifts 
Q . Large scale structure is detected in three dimensions 
without the detection of individual galaxies. This process 
is referred to as 21 cm intensity mapping. These experi- 
ments are sensitive to structures at a redshift range that 
is observationally difficult to observe for ground-based 
optical experiments due to a lack of spectral lines. Yet 
these experiments are extremely economical [6] since they 
only require limited resolution and no moving parts. 

Intensity mapping is sensitive to both the Baryon 
Acoustic Oscillations (BAO) and to weak lensing, two 
of the most powerful observational methods to deter- 
mine cosmological parameters. It has been shown that 
BAO detections from 21cm intensity mapping are pow- 
erful probes of dark energy, comparing favourably with 
Dark Energy Task Force Stage IV projects within the 
figure of merit framework 0, [1| . 

In this paper we present projected constraints on mod- 
ified gravity models from 21cm intensity mapping. In 
Section [TT] we describe the modified gravity models con- 
sidered. In Section IIIII we discuss the observational sig- 
natures accessible to 21 cm intensity mapping, and cal- 
culate the effects of modified gravity on these signatures. 
In Section HVl we present statistical analysis and results 
and we conclude in Section W\ 
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We assume a fiducial ACDM cosmology with WMAP5 
cosmological parameters: Vt m = 0.258, f2f, = 0.0441, 
fl A = 0.742, h = 0.719, n s = 0.963 and log 10 A s = -8.65. 
We will follow the convention that lo x = h 2 Vl x . 
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II. MODIFIED GRAVITY MODELS 

Here we describe some popular modified gravity mod- 
els for which projected constraints will later be derived. 
Throughout we will use units in which G = c = h = 1 
and will be using a metric with mostly negative signature: 

(+,---)• 



A. f(R) Models 

In the f(R) paradigm, modifications to gravity are in- 
troduced by changing the standard Einstein-Hilbert ac- 
tion, which is linear in R, the Ricci scalar. The mod- 
ifications are made by adding an additional non linear 
function of R (§413 



R + f(R) 

I6tt 



r 



(1) 



where L m is the matter Lagrangian. See [12| for a com- 
prehensive review of f(R) theories of gravity. 

The choice of the function f{R) is arbitrary, but in 
practice it is highly constrained by precise solar system 
and cosmological constraints, as well as stability criteria 
[ill, [HI (see below). In this paper, we choose parameteri- 
zations of f(R) such that it asymptotes to a constant for 
a certain choice of parameters and thus approaches the 
fiducial ACDM. 

In general, f(R) models have enough freedom to mimic- 
exactly the ACDM expansion history and yet still impose 
a significant modification to gravity [l5|, [l6j]. As such 
probes of the expansion history are less constraining than 
probes of structure growth, which will be evident in the 
constraints presented in later sections. 

Variation of the above action yields the modified Ein- 
stein Equations 

G^+SrR^- - Df^j g„ v -V„V„f R = 8irT^ (2) 

where = df(R)/dR, a convention that will be used 
throughout. f(R) gravity is equivalent to a scalar-tensor 
theory [TH, [T3| with the scalar field fn having a mass and 
potential determined by the functional form of f(R). The 
field has a Compton wavelength given by its inverse mass 
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The main criterion for stability of the f(R) model is that 
the mass squared of the /# field is positive, i.e. fun > 0. 
In most cases, this simply corresponds to a sign choice 
for the field fx (specifically for the model we consider 
below, fno is constrained to be less than 0). 

On scales smaller than Ac , gravitational forces are en- 
hanced by 4/3, while they reduce to unmodified gravity 
on larger scales. The reach of the modified forces Ac 



generically leads to a scale-dependent growth in f(R) 
models. 

While the dynamics are significantly changed in f(R), 
the relation between matter and the lensing potential is 
unchanged up to a rescaling of the gravitational constant 
by the linear contribution in /. The fractional change 
is of order the background field value Jr = f R (R) <C 1 
where R is the background curvature scalar. 

Proceeding further requires a choice of the functional 
form for /. A functional form is considered which is 
representative of many other cases. 

Hn and Sawicki proposed a simple functional form for 
f(R) [18], which can be written as 



f(R) = -Ro 



ci(R/R y 

c 2 (R/R ) n - 



1 



(4) 



where we have used the value of the scalar curvature 
in the background today, Rq = R\ z —o for convenience. 
This three parameter model passes all stability criteria 
for positive n, c\ and C2- One parameter can be fixed 
by demanding the expansion history to be close (within 
observational limits) to ACDM. In this case, Equation 
H] can be conveniently reparametrized and approximated 
by 



/(J2) w-2A- 



fnaRo ( -Ro 
R 



(5) 



Here A and /ro — the value of the /r field in the back- 
ground today — have been used to parameterize the func- 
tion in lieu of c\ and C2- This approximation is valid as 
long as l/ftol *C 1, which is necessary to satisfy current 
observational constraints [HI, [l|| . While A is conceptu- 
ally different than vacuum energy, it is mathematically 
identical and will thus be absorbed into the right hand 
side of the Friedmann equation and parameterized by 
f^A- In quoting constraints, we will marginalize over this 
parameter as it is of no use in identifying signatures of 
modified gravity. The parameter fno can be though of as 
controlling the strength of modifications to gravity today, 
while higher n pushes these modifications to later times. 
The effects of changing these parameters are discussed in 
greater detail in [18J. 

Allowed f(R) models exhibit the so-called chameleon 
mechanism: the /# field becomes very massive in dense 
environments and effectively decouples from matter. 
This effect is active whenever the Newtonian potential is 
of order the background fn field. Since cosmological po- 
tential wells are typically of order 10 -5 for massive halos, 
the chameleon effect becomes important if |/#| < 10~ 5 . 
If the background field is ~ 10~ 7 or smaller, a large 
fraction of the collapsed structures in the universe are 
chameleon-screened, so that the model becomes observa- 
tionally indistinguishable from ACDM. 

Since the chameleon effect will affect the formation of 
structure, standard fitting formulas based on ordinary 
GR simulations, such as those mapping the linear to the 
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nonlinear power spectrum, cannot be used for these mod- 
els. Recently, however, self-consistent N-body simula- 
tions of f(R) gravity have been_performed which include 
the chameleon mechanism |20l - |22j . We will use the simu- 
lation results for forecasts of weak lensing in the nonlinear 
regime below. 

It should be noted that f(R) models are not without 
difficulties. In particular, an open issue is the problem of 
potential unprotected singularities 23 25]. 

B. DGP Braneworld 

A theory of gravity proposed by Dvali, Gabadadze and 
Porrati (DGP) assumes that our four dimensional uni- 
verse sits on a brane in five dimensional Minkowski space 
26]. On small scales gravity is four dimensional but, 
on larger scales it becomes fully five dimensional. Here 
we parameterize DGP by r c , the scale at which gravity 
crosses over in dimensionality. The DGP model has two 
branches depending on the embedding of the brane in 
5D space. In the self-accelerating branch, the universe 
accelerates without need for a cosmological constant if 
r c ~ 1/ Hq [13, . In this branch, assuming a spatially 
flat Universe for now, the modified Friedmann equation 
is given by 



which clearly differs from ACDM. Thus, in contrast to 
the other models considered here, DGP without a cos- 
mological constant does not reduce to ACDM and it is 
possible to completely rule out this scenario (where the 
others can only be constrained). In fact DGP (without a 
cosmological constant) has been shown to be in conflict 
with current data [29] . It is presented here largely for 
illustrative purposes. 

On scales much smaller than r c , gravity is four- 
dimensional but not GR. On these scales, DGP can be de- 
scribed as an effective scalar-tensor theory [30-32]. The 
massless scalar field, the brane-bending mode, is repul- 
sive in the self-accelerating branch of DGP. Hence, struc- 
ture formation is slowed in DGP when compared to an 
effective smooth Dark Energy model with the same ex- 
pansion history. While the growth of structure is thus 
modified in DGP even on scales much smaller than r c , 
gravitational lensing is unchanged. In other words, the 
relation between matter overdensities and the lensing po- 
tential is the same as in GR [33j . 

As in f(R), the DGP model contains a nonlinear mech- 
anism to restore GR locally. This Vainshtein mechanism 
is due to self-interactions of the scalar brane-bending 
mode which generally become important as soon as the 
density field becomes of order unity. In the Vainshtein 
regime, second derivatives of the field saturate, and thus 
modified gravity effects are highly suppressed in high- 
density regions [3l|, HH, [HJ • We will only consider linear 
predictions for the DGP model here. 



III. OBSERVATIONAL SIGNATURES 

In this Section we describe the observational signatures 
available to 21 cm intensity mapping. We also give de- 
tails on calculating the observables within modified grav- 
ity models. We consider two types of measurements: 
the Baryon Acoustic Oscillations and weak gravitational 
lensing. 

For the fiducial survey, we assume a 200 m x 200 m 
cylindrical telescope, as in We will also present lim- 
ited results for a 100 m x 100 m cylindrical telescope to 
illustrate effects of reduced resolution and collecting area 
on the results. This latter case is representative of first 
generation projects @. In the 200 m case we assume 4000 
receivers, and in the 100m case 1000 receivers. We as- 
sume either telescope covers 15000 sq. deg. over 4 years. 
We assume neutral hydrogen fraction and the bias re- 
main constant with fini = 0.0005 today and 6 = 1. The 
object number density is assumed to be n — 0.03 per cu- 
bic ft, _1 Mpc (effectively no shot-noise, as should be the 
case in practice 0])- 



A. Baryonic acoustic oscillation expansion history 

test 

Acoustic oscillations in the primordial photon-baryon 
plasma have ubiquitously left a distinctive imprint in the 
distribution of matter in the universe today. This pro- 
cess is understood from first principles and gives a clean 
length scale in the universe's large scale structure, largely 
free of systematic uncertainties and calibrations. This 
can be used to measure the global cosmological expan- 
sion history through the angular diameter distance, cLa, 
and Hubble parameter, H, vs redshift relation. The de- 
tailed expansion and acceleration will differ between pure 
cosmological constant and modified gravity models. 

We use essentially the method of [Hj] for estimating 
distance errors obtainable from a BAO measurement, in- 
cluding 50% reconstruction of nonlinear degradation of 
the BAO feature. We assume the frequency range corre- 
sponding to z < 2.5 is covered (the lower z end should 
be covered by equivalent galaxy redshift surveys if not a 
21cm survey). For the sky area and redshift range sur- 
veyed, the 200m telescope is nearly equivalent to a perfect 
BAO measurement. The limited resolution and collect- 
ing area of the 100m telescope substantially degrades the 
measurement at the high-z end. 

The expansion history for modified gravity models can 
be calculated in an analogous way to that in General 
Relativity. The Friedmann Equation in DGP, Equation[6] 
can be written as 




where k is the curvature, and r c is the crossover scale. 
It is convenient to introduce the parameter w rc = l/4r^ 
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which stands in for r c . This equation can be solved nu- 
merically to calculate the observable quantities. 

We now calculate the expansion history in the HS 
f(R) model using a perturbative framework which is well 
suited for calculating constraints on /ro- Working in the 
conformal gauge and mostly negative signature, we start 
with the modified Einstein's Equation @. At zeroth or- 
der the left hand side of the 00 component contains the 
modified Friedmann equation 



H 2 = ^ + fR g n (a, a, ii, a) 



(8) 



where p is the average density (including contributions 
from A), the over-dot represents a conformal derivative 
and 



-1 



(/"+2A)a 2 



G 

'a a 3aa* 



For verifiability we quote 

a 2 Rl(2ad 2 - 7aa 2 + 2a* da) 



9i 



36a 3 



(9) 



(10) 



Evaluating Equation [8] at the present epoch yields the 
modified version of the standard constraint 



(11) 



Note that the modified version of the Friedmann Equa- 
tion is third order instead of first order, however, it has 
been shown that the expansion history stably approaches 
that of ACDM for vanishing /ro [18( . For observationally 
allowed cosmologies fno < 1 we expand 



fm9n = fmg n {a,a,a, a) +O{f R0 ) 



(12) 



where a is the solution to the standard GR Friedmann 
equation. 

By using Equation [T^] in Equation [5] and keeping only 
terms linear in /ro, the expansion history can be calcu- 
lated in the regular way, along with the observable quan- 
tities (Ia{z) and H(z). For small fno this agrees well with 
the calculation in [l8j where the full third order differen- 
tial equation was integrated 

In calculating the Fisher Matrix, this treatment is ex- 
act because the Fisher Matrix depends only on the first 
derivative of the observables with respect to the model 
parameters, evaluated at the fiducial model. 



B. Weak Lensing 

A second class of observables measures the spatial per- 
turbations in the gravitational metric. Modified gravity 
will change the strength of gravity on large scales and 
thus modify the growth of cosmological structure. Weak 



gravitational lensing, the gravitational bending of source 
light by intervening matter, is a probe of this effect. 

Weak lensing measures the distortion of background 
structures as their light propagates to us. Here, the back- 
ground structure is the 21 cm emission from unresolved 
sources. While light rays are deflected by gravitational 
forces, this deflection is not directly observable, since we 
don't know the primary unlensed 21cm sky. However, 
weak lensing will induce correlations in the measured 
21 cm background, since neighbouring rays pass through 
common lens planes. While the deflection angles them- 
selves are small (of order arcseconds) the deflections are 
coherent over scales of arcminutes. In this way, the lens- 
ing signal can be extracted statistically using quadratic 
estimators 36]. Given the smallness of the lensing ef- 
fect, a high resolution (high equivalent number density 
of "sources") is necessary to detect the effect. 

The weak lensing observable that is predicted by the- 
ory is the power spectrum of the convergence k. It is 
given by 



X Xa z (x) 



(13) 

where x denotes comoving distances, P{k, x) is the (lin- 
ear or nonlinear) matter power spectrum at the given 
redshift, and we have assumed flat space. The lensing 
weight function Wl(x) i s given by: 



W L ( X ) 



*(X) 



y dN 
cb^W^-X) _(,.)• d4) 



Here, dN/dz is the redshift distribution of source galax- 
ies, normalized to unity. The factor e(x) i n Equation 1131 
encodes possible modifications to the Poisson equation 
relating the lensing potential to matter (Section [H]). In 
f(R), it is given by e{x) = (1 + /fi(x))" 1 , while e = 1 
for GR as well as DGP. Note that for viable f(R) mod- 
els, e — 1 < 0.01, so the effect of e on the lensing power 
spectra is very small. 

The CAMB Sources module [13, [H[ was used to cal- 
culate the lensing convergence power spectrum in flat 
ACDM models. The HALOFIT [H interface for CAMB 
was used for calculations that include lensing at nonlinear 
scales. 

For the modified gravity models in the linear regime, 
the convergence power spectra were calculated using the 
Parametrized Post-Friedmann (PPF) approach [40j] as 
in Briefly, the PPF approach uses an interpolation 
between super-horizon scales and the quasi-static limit. 
On super- horizon scales (k <C aH), specifying the back- 
ground expansion history, together with a relation be- 
tween the two metric potentials, already determines the 
evolution of metric and density perturbations. On small 
scales (k ^> aH), time derivatives in the equations for 
the metric perturbations can be neglected with respect 
to spatial derivatives, leading to a modified Poisson equa- 
tion for the metric potentials. The PPF approach uses a 
simple interpolation scheme between these limits, with a 
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few fitting parameters adjusted to match the full calcula- 
tions [40| ■ The full calculations are reproduced to within 
a few percent accuracy. We use the transfer function of 
(4lj to calculate the ACDM power spectrum at an ini- 
tial redshift of Zi = 40, were modified gravity effects are 
negligible, and evolve forward using the PPF equations. 

For the f(R) model, we also calculate predictions in 
the nonlinear regime. For these, we use simulations of 
the HS model with n — 1 and fno values ranging from 
1(T 6 to 1CT 4 . We use the deviation AP(k)/P(k) of the 
nonlinear matter power spectrum measured in f(R) sim- 
ulations from that of ACDM simulations with the same 
initial conditions plj . This deviation is measured more 
precisely than P(k) itself. We then spline-interpolate the 
measurements of AP(k)/P(k) for k = 0.04 - 3.1 h/Mpc 
and at scale factors a = 0.1, 0.2,. ..1.0, and multiply the 
standard nonlinear ACDM prediction (HALOFIT) with 
this value. For values of k > 3.1 h/Mpc, we simply set 
AP(k) — 0. However, for the angular scales and redshifts 
considered here (£ < 600, see below), such high values of 
k do not contribute significantly. 

One might be concerned that this mixing of methods, 
for calculating the lensing spectrum, might artificially ex- 
aggerate the effects of modified gravity if these methods 
do not agree perfectly. While the spectra calculated for 
the fiducial ACDM model differed by up to a percent 
between these methods, presumably due to slight differ- 
ences in the transfer function, this should have no effect 
on the results. Any direct comparison between spectra 
(for example finite difference derivatives) are made be- 
tween spectra calculated in the same manner. Note that 
the Fisher Matrix depends only on the first derivative of 
the observables with respect to the parameters and no 
cross derivatives are needed. 

The lensing spectra were not calculated for non-flat 
models, but it is expected that the CMB and BAO are 
much more sensitive to the curvature and as such the 
lensing spectra are relatively unaffected. Formally we 
are assuming that 



Reconstructing weak lensing from 21 cm intensity maps 
involves the use of quadratic estimators to estimate the 
convergence and shear fields. The accuracy with which 
this can be done increases with information in the source 
maps, however, this information saturates at small scales 
due to nonlinear evolution. As such, one cannot improve 
the lensing measurement indefinitely by increasing reso- 
lution, and the experiments considered here extract much 
of the available information within the redshift range con- 
sidered. 

The accuracy with which the convergence power spec- 
trum can be reconstructed from 21 cm intensity maps was 
derived in [3|| , where the effective lensing galaxy density 
was calculated at redshifts 1.25, 3 and 5 (see Figure 7 and 
Table 2 therein) . The effective volume galaxy density was 
corrected for the finite resolution of the experiment con- 



sidered here. It was then interpolated, using a piecewise 
power law, and integrated from redshift 1 to 2.5 to obtain 
an effective area galaxy density of n g /o-\ = 0.37arcmin -2 . 
The parameter a\ is the variance in the intrinsic galaxy 
ellipticity, which is only used here for comparison with 
optical weak lensing surveys. From the effective galaxy 
density the error on the convergence power is given by 



A ° KK W = \ op+^f ( C ' KK W + — ) ( 15 ) 

V (2e+l)fsky V Tig J 

where f s k y is the fraction of the sky surveyed. The galaxy 
distribution function dN/dz used to calculate the theo- 
retical curves (from Equation IT3|) should follow the ef- 
fective galaxy density. Instead for simplicity, a flat step 
function was used, with this distribution function equal 
from redshift 1 to 2.5 and zero elsewhere. While the dif- 
ference between the these distributions would have an 
effect on the lensing spectra, the effect on differences of 
spectra when varying parameters is expected to be neg- 
ligible. Our approximation is also conservative, since the 
proper distribution function is more heavily weighted to- 
ward high redshift. Rays travelling from high redshift 
will be affected by more intervening matter and thus ex- 
perience more lensing. This would increase the lensing 
signal, allowing a more precise measurement. 

Figure Q] shows the lensing spectra for the fiducial cos- 
mology and a modified gravity model, including both lin- 
ear and nonlinear calculations. The linear regime is taken 
to be up to I = 140 for projected constraints. For cal- 
culations including weak lensing in the nonlinear regime, 
C KK (£) up to I — 600 is used for the larger telescope. 
Beyond this scale the model used for lensing error-bars 
is not considered accurate at the shallowest redshifts in 
the source window [36} . This cut off coincides with the 
scale at which information in the source structures sat- 
urates due to non-linear evolution in standard gravity 
(although it is also not far from the resolution limit of 
the experiment). We speculate that a similar phenomena 
would occur in modified gravity and smaller scales are 
not expected to carry significant additional information. 
Note that it is the source structures in which information 
saturates. At smaller scales the lensing spectrum would 
continue to carry information [42j if it could be recon- 
structed. For the smaller telescope the scale is limited 
to £ < 425 by the resolution at the high end of the red- 
shift window. If the redshift window were subdivided into 
narrower bins, it would be possible to use information at 
scales down to I f=s 1000 in the centre bins as at these red- 
shifts the telescope resolutions are better and structures 
are less non-linear. However, considering tomographic 
information is beyond the scope of this work. It is noted 
that these scales are very large by weak lensing standards 
where optical surveys typically make detections down to 
an I of order 10 5 . 
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FIG. 1: The Weak lensing convergence power spectra for 
ACDM and the HS f(R) model with n = 1 and f R0 = 1(T 4 . 
Galaxy distribution function is flat between z = 1 and z = 
2.5. 

C. External Priors from Planck 

While the CMB is not sensitive to the late time effects 
of modified gravity (except by the integrated Sachs- Wolfe 
effect), it is invaluable for constraining other parameters 
and breaking degeneracies. As such, projected informa- 
tion from the Planck experiment is included. The Planck 
covariance matrix used here is given in [i^, Table II] . All 
late time cosmological parameters (including the curva- 
ture) are marginalized over, removing information con- 
tained in the ISW effect, and ensuring that sensitivity to 
f(R) is entirely from 21 cm tests below. The only remain- 
ing parameter that is related to the late time expansion 
is 9 S , the angular size of the sound horizon, which is then 
used as a constraint on the parameter sets of the modified 
gravity models. 

IV. RESULTS 

To quantify the projected constraints on f(R) models, 
the Fisher matrix formalism is employed. The HS f(R) 
models reduces to the fiducial model for vanishing ffto 
and any value of n. Thus the Fisher Matrix formalism 
is used to project constraints on for given values of 
n. In the case of DGP, which does not reduce to the 
fiducial model, it is shown that a measurement consis- 
tent with the fiducial model can not be consistent with 
DGP for any parameter set. Unless otherwise noted, we 
account for freedom in the full cosmological parameter 
set: h, co m , w&, ^>k, A s and n s ; representing the Hubble 
parameter; physical matter, baryon and curvature densi- 
ties; amplitude of primordial scalar fluctuations and the 
spectral index; respectively. 

Within the f(R) models, the fiducial model is a special 
point in the parameter space as there are no modifica- 



tions to gravity. As such, one cannot in general expect 
perturbations to observables to be linear in the f(R) pa- 
rameter fno, an assumption implicit in the Fisher Ma- 
trix formalism. This assumption does seem to hold for 
the expansion history, where our first order perturbativc 
calculation agrees with the full solution to the modified 
Friedmann Equations calculated in [l8j]. However, this is 
not the case for weak lensing. For each f(R) model, the 
lensing spectrum was calculated for several values of fm. 
It was observed that enhancements to the lensing power 
spectrum go as 

with a(£) in the 0.5-0.7 range. This is because the reach 
of the enhanced forces in f(R) is a power law in Jrq fol- 
lowing Equation ([3]), and the enhancement of the power 
spectrum for a given mode k roughly scales with the time 
that this mode has been within the reach of the enhanced 
force. Because of this behaviour, the constraints derived 
within the Fisher Matrix formalism depend on the step 
size in fn used for finite differences. 

To correct for this, we use a step size that is dependent 
on the final constraint. The weak lensing Fisher Ma- 
trices where calculated for fn step sizes of 10~ 3 , 10~ 4 
and 10 -5 . These were then interpolated — using a power 
law — such that the ultimate step size used for finite dif- 
ferences is roughly the quoted constraint on the modified 
gravity parameter. For instance when the 95% confidence 
constraint on //jo is quoted, the step size for finite differ- 
ences is A//jo ~ 2<7 f R0 , where <jj ro is calculated from the 
interpolated Fisher matrix. This is expected to be valid 
down to step sizes at the 10~ 6 level where the chameleon 
mechanism is important. As such, for constraints below 
10 -5 a step size of 10 -5 is always used. Note that this is 
conservative because an over sized finite difference step 
always underestimates the derivative of a power law with 
an power less than unity. For constraints above the 10 -3 
level a step size of 10~ 3 is used, which is the largest mod- 
ification to gravity simulated. These constraints are con- 
sidered unreliable due to these difficulties. We reiterate 
this this only affects results that include weak lensing in- 
formation. Likelihood contours remain perfect ellipses in 
this procedure (which is clearly inaccurate), however the 
spacing between contours at different confidence levels is 
altered. 

Figure [5] shows the projected constraints on the HS 
f(R) model with n = 1 for various combinations of ob- 
servational techniques, and a (200m) 2 telescope. The 
elements in the lensing fisher matrix associated with the 
curvature are taken to be zero for the reasons given in 
Section UlI Bl While this assumption is not conservative, 
it is expected to be valid, as the angular diameter dis- 
tance as measured by the BAO is very sensitive to the 
curvature. In total three f(R) models were considered: 
HS with n — 1,2,4. The results are summarized in Table 
III 

It was found that while weak lensing, in the linear 
regime, is very sensitive to the modifications to grav- 
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FIG. 2: Projected constraints on the HS f(R) model with 
n = 1 using several combinations of observational techniques, 
for a 200 m telescope. All curves include forecasts for Planck. 
Allowed parameter values are shown in the /ho — h plane at 
the 68.3%, and 95.4% confidence level. Results are not shown 
for "WL" which were calculated much less accurately (see 
text). 



95% confidence 
upper limits 


HS 


fm\ 


n = 1 


NL WL 


n = 2 


n = 4 


BAO 
WL 

BAO+WL 


1.5e-02 
2.3e-03 
5.0e-05 


4.3e-05 
8.9e-06 


1.8e-02 
4.0e-03 
9.7e-05 


3.0e-02 
8.6e-03 
4.6e-04 



TABLE I: Projected constraints on f(R) models for various 
combinations of observational techniques, for a 200 m tele- 
scope. Constraints are the 95% confidence level upper limits 
and include forecasts for Planck. The non linear results (col- 
umn marked NL WL) are for the HS model with n = 1. 
Results that make use of weak lensing with constraints above 
10 -3 are only order of magnitude accurate. The linear regime 
is taken to be I < 140, with the nonlinear constraints extend- 
ing up to I = 600. 



ity, it is only barely capable of constraining f(R) models 
without separate information about the expansion his- 
tory. Even with the inclusion of Planck forecasts, degen- 
eracies with h and the mean curvature, drastically 
increase the uncertainties on the modified gravity param- 
eters. Indeed these three parameters are more than 95% 
correlated (depending on the exact model and confidence 
interval) . This of course brings into question the neglect 
of the uJk terms in the weak lensing Fisher Matrix. How- 
ever it is noted that in these cases, the predicted limits 
on the curvature are \u>k\ < 0.025 at 95% confidence. 
The current, model independent, limits on the curvature 
using WMAP, SDSS and HST data are approximately 
half this value (44]. Our neglect of any direct probes 
of the expansion history for the Planck+WL constraints 
is clearly unrealistic; however, the constraints illustrate 
what is actually measured by weak lensing. In any case 



FIG. 3: Same as Figure[2]but for a 100 m cylindrical telescope. 

these degeneracies are broken once BAO measurements 
are included, and in this final case the modified gravity 
parameters are correlated with the other parameters by 
at most 35%. Also, considering lensing in the nonlinear 
regime breaks the degeneracy to a certain extent. 

First generation cylindrical telescopes will likely be 
smaller than the one considered above. To illustrate the 
differences in constraining ability, we now present a few 
results for a cylindrical radio telescope that is 100m on 
the side. Reducing the resolution of the experiment de- 
grades measurements in a number of ways. BAO mea- 
surements become less than ideal in the higher redshift 
bins. The smallest scale that can be considered for weak 
lensing drops to about I = 425. A more important effect 
is that the lensing spectra can not be as accurately re- 
constructed, dropping the effective galaxy density down 
to n g ja\ = 0.22. Figure [3] shows analogous results to [2] 
but for a telescope with half the resolution. 

To show that a set of measurements consistent with the 
fiducial model would be inconsistent with DGP we first 
fit DGP to the fiducial model's CMB and BAO expansion 
history by minimizing 

X 2 = (r DGP - r flducial ) T C~ 1 (r DGP - r f lducial ) (16) 

where r^GP and r fiducial are vectors of observable quan- 
tities as calculated in the DGP and fiducial models, and 
C is the covariance matrix, r includes BAO cIa(z) and 
H(z) as well as Planck priors on uj m and 8 S . Note that 
X 2 is not truly chi-squared since r fiducial contains fidu- 
cial model predictions and is not randomly distributed 
like a real data set. 

Performing the fit yields DGP parameters: h — 0.677, 
u m = 0.112, u k = -0.0086 and uj rc = 0.067. Figured 
shows the deviation of H and cIa respectively for the best 
fit DGP model compared to the fiducial model. \ 2 — 
332.8 for the fit despite there only being 16 degrees of 
freedom, and as such a measurement consistent with the 
fiducial model would thoroughly rule out DGP. 
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FIG. 4: Ratio of the coordinate cIa(z) (top) and the Hubble 
parameter H(z) (bottom) as predicted by the best fit DGP 
model to the fiducial model. Error bars are from 21 cm BAO 
predictions. Fit includes BAO data available from the 200 m 
telescope and CMB priors on 8 3 and w m . 
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FIG. 5: Weak lensing spectra in for DGP and a smooth dark 
energy model with the same expansion history. DGP param- 
eters are h = 0.665, uj m = 0.116, uik = and uj tc = 0.06. 
Errorbars represent expected accuracy of the 200 m telescope. 



In the case that expansion history measurements are 
consistent with DGP, the question arises as to whether 
DGP could be distinguished from a smooth dark energy 
model that had the same expansion history. The addi- 
tional information in linear perturbations as measured by 
weak lensing allows DGP to be distinguished even from 
a dark energy model with an identical expansion history. 
Figure [5] shows the lensing spectra for a DGP cosmology 
similar to the best fit discussed above, as well as the dark 
energy model with the same expansion history as in [2!| . 

In principle one should consider the small amount of 
freedom within the DGP parameter set that could be 



used to make the DGP spectrum better fit the dark 
energy spectra. However this is unlikely to signifi- 
cantly change the spectrum as all relevant parameters 
are tightly constrained by the CMB and BAO. For ex- 
ample it is clear from Figure [5] that the lensing spectra of 
the two models would better agree if the amplitude of pri- 
mordial scaler perturbations was increased in the DGP 
model. However, Planck measurements would only allow 
of order half a percent increase while the disagreement 
is of order 10%. This is somewhat justified by the lack 
of correlations found in the f(R) fisher matrices once all 
three observational techniques are included. In addition 
we have not considered information from weak lensing 
in the nonlinear regime. Adding nonlinear scales would 
only make our conclusion that DGP and smooth dark 
energy are distinguishable with these observations more 
robust. 



V. DISCUSSION 

We have shown that the first generation of 21 cm inten- 
sity mapping instruments will be capable of constraining 
the HS f(R) model (with n = 1) down to a field value of 
\fno\ < 2 x 10" 5 at 95% confidence (Figure [3). This is 
an order of magnitude tighter than constraints currently 
available from galaxy cluster abundance [l9f. Further- 
more, model parameters in this regime are not ruled out 
by Solar System tests. 

In comparing Figures [2] and [3] it is clear that a more 
advanced experiment, with resolution improved by a fac- 
tor of two, would further half the allowed value of \f mi- 
lt should be noted however, that halving of the allowed 
parameter space does not correspond to a factor of four 
increase in information. Deviations in the lensing spec- 
trum scale sub- linearly in the f(R) parameters, enhanc- 
ing the narrowing of constraints as information is added 
(see Section [TVf . 

While we have concentrated on a particular f(R) 
model, many viable functional forms for f(R) have been 
proposed in the literature [45H47| . The predictions for 
the growth of structure in these different models agree 
qualitatively: the gravitational force is enhanced by 4/3 
within Ac, enhancing the growth on small scales. How- 
ever, there are quantitative differences in the model pre- 
dictions due to the different evolution of Ac over cosmic 
time. Our results for the HS model with different val- 
ues of n should thus cover a range of different functional 
forms for f(R). Table U shows that our constraints do 
not depend very sensitively on the value of n. This is be- 
cause the weak lensing measurements cover a wide range 
of scales as well as redshifts. Furthermore, it is straight- 
forward to map the enhancement in the linear P(k, z) at 
given k and z from the HS model considered here to any 
other given model, to obtain approximate constraints for 
that model. 

Future cluster constraints will almost certainly im- 
prove on the current limits of |/ro| < few 10~ 4 [lj|. 
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However, for smaller field values, the main effect of f(R) 
gravity shifts to lower mass halos, since the highest mass 
halos are chameleon-screened (see Fig. 2 in [22]). Hence, 
future cluster constraints will depend on the ability to ac- 
curately measure the halo abundance at masses around 
few 1O 14 M and less. Further more, the constraints from 
cluster abundances depend sensitively on the knowledge 
of the cluster mass scale, and are already systematics- 
dominated [jj}. Weak lensing constraints have a com- 
pletely independent set of observational systematics, and 
are in principle less sensitive to baryonic or astrophysi- 
cal effects. Thus, the forecasted constraints on modified 
gravity presented here are quite complementary to con- 
straints from cluster abundances. 

The processes that produce the BAO feature in the 
matter power spectrum are understood from first princi- 
ples. In addition the BAO length scale can be extracted 
even in the presence of large uncertainties in biases and 
mass calibrations. Likewise, weak lensing on large scales 
is well understood, with baryonic physics being much 
less important than on smaller scales [Hj]. In addition 
the dominant systematics present in optical weak lens- 
ing surveys are instrumental in nature and not intrinsic 
to the quantities being measured. While 21 cm intensity 
mapping is as yet untested, instrumental systematics will 
be very different from those that affect the optical. 

In the case of this study, and more generally for cosmo- 
logical models which substantially modify structure for- 
mation, the motivation for higher resolution comes not 
from improved BAO measurements but from better weak 
lensing reconstruction. Higher resolution not only makes 
weak lensing information available at higher multi-poles, 
but improves the accuracy at which lensing can be recon- 
structed on all scales. 

The inclusion of lensing information in the nonlinear 
regime was crucial, and largely responsible for the com- 
petitiveness of these forecasts. As seen in Figure [TJ much 
of the constraints come from multi-poles in the nonlinear 
regime. It should be noted that for the higher resolution 
experiment considered, the minimum scale is limited not 
by the resolution at high redshift, but by the saturation 
of information in nonlinear source structures at low red- 
shift m. 

Our constraints from lensing are conservative since 
only one wide source redshift bin was considered, limited 
to £ < 600 as described above. To maximize informa- 
tion, the source redshift range could be split into multiple 
bins, properly considering the correlation in the lensing 
signal between them; a process known as lensing tomog- 
raphy. The low redshift bin would be limited as above, 
and the high redshift bin would be limited by the resolu- 
tion to £ ~ 850 at z = 2.5. However in intermediate bins, 
the lensing signal could be reliably reconstructed above 
<?« 1000. 



Unlike most smooth dark energy models, such as 
quintessence, constraints on the models considered here 
are chiefly sensitive to structure formation, as is clear 
from Figure [5] No experiment in the foreseeable fu- 
ture will be able to improve upon the constraints we 
project on these models using exclusively expansion his- 
tory probes (except by breaking degeneracies). These 
forecasts show that 21 cm intensity mapping is not only 
sensitive to a cosmology's expansion history through the 
BAO, but also to structure growth through weak lensing. 
The weak lensing measurements cannot compete with far 
off space based surveys like Euclid or JDEM, which will 
have galaxy densities of order 100arcmin _2 Q and reso- 
lution to far greater I. However, cylindrical 21 cm exper- 
iments are realizable on a much shorter time scale and 
at a fraction of the cost. In addition, the measurements 
considered here are approaching the limit at which f(R) 
models can be tested. For |/ro| much less than 10^ 6 the 
chameleon mechanism becomes important before there 
are observable modifications to structure growth, reduc- 
ing the motivation to further study these models. 

It has also been shown that, for these experiments, a 
BAO measurement consistent with ACDM would defini- 
tively rule out DGP without a cosmological constant as 
a cosmological model. Even in the case that a BAO mea- 
surement consistent with DGP is made, the model is still 
distinguishable from an exotic smooth dark energy model 
through structure growth. The former result is not sur- 
prisin g g iven that DGP is now in conflict with current 
data [23. However it is illustrative that a single exper- 
iment can precisely probe both structure formation and 
expansion history. Even a dark energy model that con- 
spires to mimic DGP is, to a large extent, distinguishable. 

We have studied the effects of modified gravity theories 
on observational quantities for future 21 cm surveys. Be- 
cause these surveys measure the distribution of galaxies 
on large angular scales over large parts of the sky, they 
are well suited to measure the expected deviations rel- 
ative to standard general relativity. We have computed 
the predictions of modified gravity in the linear and non- 
linear regimes, and compared to the sensitivity of future 
surveys. We find that a large part of parameter space 
can be tested. 
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